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Abstract. Ordering by thermal fluctuations is studied for the classical XY 
antiferromagnet on a checkerboard lattice in zero and finite magnetic fields by 
means of analytical and Monte Carlo methods. The model exhibits a variety of 
novel broken symmetries including states with nematic ordering in zero field and 
with triatic order parameter at high fields. 
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1. Introduction 

A huge degeneracy of classical ground states in geometrically frustrated magnets can 
be lifted by quantum or thermal fluctuations via a so called order by disorder effect 
PQ. In the present work we study the thermal order by disorder effect for the XY 
antiferromagnet on a checkerboard lattice. This lattice is a two-dimensional network of 
corner-sharing squares with crossings, which are topologically equivalent to tetrahedra, 
see Fig. ^ The present model can, therefore, be relevant to real pyrochlores with the 
easy-plane type anisotropy Ei^T^Oy and Ei'2Sn207 [5] and also have an experimental 
realization as an array of the Josephson junctions or a superconducting wires network 
in transverse magnetic field 

The considered model is described by the Hamiltonian 



where Si = (cos ifi, sin. <pi) are the classical planar spins and H = (H,0) is in- plane 
magnetic field. The sum is over nearest-neighbors pairs on a checkerboard lattice and 
J > is an antiferromagnetic coupling constant. The Hamiltonian can be expressed 
as a sum over elementary plaquettes (squares with crossings): 



where S p = ^2 i£p Si is the total magnetisation of a plaquette, N is the number of 
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The minimal total energy is reached if the above constraint is satisfied on every 
plaquette. In the field range < H < H sat = 8 J the ground state of the model remains 
underconstrained and infinitely degenerate with a finite entropy. This huge degeneracy 
is determined by two factors: (i) possible basis spin quartets, which obey the classical 
constraint, are parametrized by two continuous variables and (ii) by various periodic 
and aperiodic spin structures formed with the same basis quartet. The latter property 
can be seen by looking at the Fourier transform of the interactions on the checkerboard 
lattice 0], which has a flat momentum-independent branch with minimal energy. 

2. Thermal order by disorder effect 

Since the degeneracy of the classical ground state is a consequence of lattice topology 
rather than being a symmetry imposed property, various classical ground state 
configurations have different excitation spectra. At finite temperatures magnetic 
system spans a phase volume in the vicinity of the ground state manifold. Due 
to a varying density of excitations, the system can be effectively trapped in the 
neighbourhood of certain ground states. Such an ergodicity breaking leads to a 
lifting of zero-T degeneracy of a frustrated magnet and to thermal order by disorder 
selection. The statistical weights of different ground state configurations ip are given 
by w[ip] ~ exp(— F[ip]/T), where F[ip] is the partial free energy obtained by integrating 
out 'fast' excitation modes. The minimum of F[ip] ensures that a macroscopic system 
is trapped in the vicinity i/> mm and probability to find it in another classical ground 
state is vanishingly small. 

To find which spin configurations are favoured by thermal fluctuations we start 
with a simple perturbative treatment over the mean-field result. In the ground state 
configuration a magnitude of a local field derived from Eq. @ is the same H\ oc = 2 J 
for all sites and all external fields < H < H sat ■ The harmonic spin- wave Hamiltonian 
is expressed as 

n 2 = -tfioc Yl s t + J Y1 s f s j cos e n > s i » 1 - \ s f > ( 4 ) 

» (ij) 

where components of every spin are taken in its local coordinate frame and 9ij is an 
angle between neighbouring spins. The first term in H.2, which describes uncorrelated 
fluctuations of individual spins with (S^ 2 ) — T/H\ OCl is taken as an unperturbed spin- 
fluctuation Hamiltonian, whereas the second term is a perturbation V. The correction 
to the free-energy is given by AF = —(V 2 )/2T: 

AF = -(T/8) J2 cos2 % - -(T/8) J2 (S, • S,) 2 . (5) 

(v) (ij) 

Thermal fluctuations produce, therefore, an effective biquadratic exchange, which lifts 
the zero-T degeneracy in favour of maximally collinear states with the largest number 
of co&Oij = ±1. The harmonic Hamiltonian can be, of course, diagonalized with the 
help of the Fourier transform to obtain the spin wave modes u£ and their contribution 
to the free-energy: 

AF 2=T^£lnKAT). (6) 

n k 

The problem is reduced to minimization of the sum in JSJ. Still the real-space 
second order perturbation result appears to give a correct first insight especially for 
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multisublattice configurations with n > 2, when exact diagonalization of TL2 becomes 
increasingly cumbersome. 

An additional complication of highly frustrated magnets stems from a presence of 
several branches of zero modes: w£ l = 0, in which case Eq. © is not anymore correct. 
Instead the low-T contribution to the free energy becomes: 

AF = (iV 2 / 2 + N 4 /4)T ln( J/T) +T ^ lnw£/J + T/ 4 , (7) 

where N2 is the number of usual harmonic or quadratic modes, A/4 is the number of 
zero or soft modes, and /j is a contribution from interaction between soft modes and 
their interaction with quadratic modes. In order to estimate the last term one has 
to solve a nonlinear problem, which is, generally, a very complicated task. However, 
a partial selection between various spin configuration can be done on a basis of the 
leading T\n{J/T) term. If the total number N 2 + A 4 is fixed (= N for the XY 
checkerboard antiferromagnet), then, the free energy is minimal for states with the 
maximum number of soft modes. The number of soft modes can be found either from 
direct diagonalization of H2 or from geometric consideration, which assigns a local soft 
mode to every void (empty square) with all spins around it being parallel or antiparallel 
to each other [SJE10- Thus, the soft modes act similar to harmonic excitations and 
stabilize collinear states. The presence of soft modes is most easily seen in the low- 
temperature specific heat, which is found from Eq. Q to be C = \ — N4/AN. The 
deviation of the specific heat from a universal value of C = | tells how many soft 
modes exist in a low-temperature state of the XY checkerboard antiferromagnet. 

3. Zero field behaviour 

At zero magnetic field the ground state constraint Q specifies configurations with 
S p = on every plaquette. Such configurations can be constructed either from 
noncollinear or collinear spin quartets. As was argued above, thermal fluctuations tend 
to select maximally collinear states with two up- and two down-spins on an arbitrarily 
chosen axis in the XY plane. The gauge transformation Sf down — — 5,f down maps H2 
for an arbitrary collinear state on the same reference harmonic Hamiltonian. Every 
collinear state has, therefore, the same harmonic spectrum and the same number of 
zero (soft) modes A4 = ~JV. Hence, the specific heat of a collinear state is C = |. All 
terms in the low-T expression for the free energy Q except the last one coincide for 
all collinear states. Our estimate indicates that the Neel state with the ordering wave- 
vector Q = (71", 7r) on the original square lattice has the lowest anharmonic contribution 
fi among all translationally symmetric states. This, however, does not necessarily 
mean an appearance of a quasi long-range order at q = Q when T — > 0. So called 
weathervane defects jS] or wondering (rough) domain walls [S] have been considered 
as a source of disorder for the kagome antiferromagnet. For the XY checkerboard 
antiferromagnet there are no zero-width domain walls and only weathervane defects 
can destroy a quasi long-range translational order. They cost zero classical energy 
and increase the free-energy by AF^ ~ e^T . A finite T-independent density of such 
defects can be estimated as na — (l + e £d ) _1 by neglecting interaction between defects. 
These defects can destroy the true long-range order if their concentration exceeds 
the percolation threshold on the corresponding lattice. In this case only nematic 
correlations desribed by a traceless second-rank tensor 0"f m = (5f S?) — ^S a/3 will be 
present at low temperatures. 
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Figure 1. Top left panel: the checkerboard lattice. Two right panels: 
temperature dependence of squares of nematic (top) and Neel (bottom) order 
parameters for three lattice sizes. Left bottom panel: finite-size scaling for the 
nematic order parameter. Temperatures are given in units of J. 



In the absence of analytical theory to deal with this sort of behaviour we have 
tried to derive the necessary information from Monte Carlo simulations. We have 
calculated squares of the two relevant low-T order parameters: 

s(Q) = ^Bsfsj> lQ(r, ~ r3) , Ci^=±Y^(S?S?S?S?)-± (8) 

The results are presented in Fig. ^ There is a clear signature of the nematic order 
seen by the enhancement of 0^ om (N) at low temperatures. The Neel order parameter 
reaches only very small values and does not show any appreciable enhancement. 
Considering a change in the scaling behaviour of the nematic order parameter 
0^ cm (N) ~ 1/N a from a = 1, short-range correlations, at high temperatures to 
a < 1, power-law correlations, at low temperatures we estimate the Kosterlits-Touless 
transition temperature as Tkt = 0.014(2) J. 



4. Finite field phases 

In external magnetic fields < H < H sat the ground states with intermediate 
magnetisation © are, generally, noncollinear, except for H = ^H sat , when a collinear 
'uuud' state, Fig. [2), belongs to the ground state manifold. Selection of the uuud states 
by thermal fluctuations leads to a 1/2-magnetization plateau, which is similar to a 
1/3-plateau of a classical kagome antiferromagnet [7]. At all other fields only a partial 
collinearity is possible in the classical ground state. Geometric consideration suggests 
two prime candidate states shown in Fig. |3 The first canted state, Fig. [5^, exists 
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Figure 2. Top left panel shows three ground state configurations with maximal 
collinearity. Top right panel presents a translational pattern with maximum 
number of soft modes for the partially collinear state (c). Bottom left panel: 
field dependence of the specific heat for N = 1024 cluster at T = 0.002 J. Bottom 
right panel: the evolution of the triatic order parameter for the same cluster. The 
inset shows the field dependence of S yy (q) corresponding to state (d) for lattice 
sizes N = 256, 400, 1024 (from top to bottom). Magnetic field is given in units 
of J. 



in the whole range of fields and does not break remaining spin reflection symmetry 
about the field direction. The second partially collinear state, Fig. with three 
identical sublattices appears only for H > ^H sat and does break the mirror symmetry. 
It is easy to check that the fluctuation induced biquadratic exchange JSJ favours a 
'more' collinear state with broken reflection symmetry AF = — (T/8)[3 + |(8/i 2 — 5) 2 ], 
h = H I H sa t over a 'less' collinear canted state AF = -(T/8)[2 + A(2h 2 - l) 2 ] in the 
whole range of existence of the former state. 

The translational degeneracy of the high-field partially collinear states is similar 
to the degeneracy of the uuud states and corresponds to the dimer coverings of a 
square lattice the total number of such states being ~ 1.157^. The problem 
of lifting translational degeneracy for the uuud states is similar to the zero-field 
case discussed above. Also, the uuud have the same specific heat C = |. For the 
partially collinear states there are additional features in the thermal order by disorder 
effect. A general partially collinear state does not have soft modes. Soft modes, 
corresponding to all parallel spins around an empty square, exist, nevertheless, for 
certain translational patterns. The maximum number of soft modes -/V4 = N/A appears 
for a state shown in Fig. EJl, which corresponds to a columnar arrangement of effective 
dimers. Such a translational pattern has the lowest 7Tn( J/T) contribution to the free 
energy. However, for all temperatures T > 0.001J accessible with our Monte Carlo 
code we did not find a nonvanishing value of the structure factor corresponding to the 
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this state. Thus, there is no a conventional Ising order parameter (Sf) = 0, Vi in the 
high-field partially collinear state. Instead an Ising reflection symmetry is broken by 
a unique triatic order parameter: 

T yyy = (SfSfS!) . (9) 

For Monte Carlo simulations of the XY chekcer-board antiferromagnet we have 
used the standard Metropolis algorithm discarding ~ 10 5 Monte Carlo steps per spin 
(MCS) to reach thermal equilibrium and, then, average observables over 10 6 -10 7 MCS. 
The lattice sizes were up to N = 1024. In order to estimate statistical errors, all 
results have been averaged over 10-20 runs. The field dependence of the specific 
heat C is shown in Fig. [21 for T — 0.002 J. C starts at value, which is very close to 
analytically predicted | . In applied field the system looses gradually soft modes up to 
H = Q.5H sat , where uuud state with C = | appears. The peaks in the specific heat 
at H c i = 4.15J and = 6.8J indicate the phase transitions to a high-field state 
with broken reflection symmetry. The order parameters for this state are presented 
in the bottom right panel of Fig. |21 The square of the triatic order parameter T yyy 
is nonzero between H c \ and H C 2 . It shows a good statistical averaging and very little 
finite size dependence: results for smaller clusters fall on top of the presented data for 
N = 1024. The inset shows S vv (q) with q = (tt,tt) on a new lattice built from one 
sort of tetrahedra, which should be nonzero if the long-range order corresponding to 
Fig. Eli is present. The obtained data show that S yy (q) = in the thermodynamic 
limit. Thus, the high-field state is described by the triatic order parameter T yyy , which 
makes the checkerboard lattice to be similar to a classical kagome antiferromagnet [7] . 
Our results suggest that the triatic state survives up to T ~ 0.015 J. 

In conclusion, we have investigated the low temperature phases of the XY 
checkerboard antiferromagnet in strong external fields. We have found that thermal 
fluctuations stabilize interesting phases with new type of broken symmetries: nematic 
order at H = and triatic order © in the field range \H sa± < H < H sat . Further 
theoretical investigations should focus on a unique H-T diagram of this model. 
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